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We consider a mass-asymmetric electron and hole bilayer. Electron and hole Coulomb correlations and electron 
and hole quantum effects are treated on first princles by path integral Monte Carlo methods. For a fixed layer 
separation we vary the mass ratio M of holes and electrons between 1 and 100 and analyze the structural 
changes in the system. While, for the chosen density, the electrons are in a nearly homogeneous state, the hole 
arrangement changes from homogeneous to localized, with increasing M which is verified for both, mesoscopic 
bilayers in a parabolic trap and for a macroscopic system. 
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1 Introduction 

Strongly correlated Coulomb systems are of growing interest in many fields, including plasmas and condensed 
matter, see e.g. [1 1 for an overview. In particular, Wigner crystal formation is one of the most prominent correla- 
tion phenomena observed in ultracold ions |2 |, dusty plasmas [3, 4|, quantum dots, e.g. fSllH and other confined 
(non-neutral) systems. Recently crystal formation in two-component (neutral) quantum plasmas was demon- 
strated by simulations |7| confirming early predictions of hole crystallization in semiconductors by Halperin and 
Rice t8J, Abrikosov [9| and others. Interstingly, this is essentially the same physical phenomenon as crystalliza- 
tion of nuclei in White Dwarf stars ifTOl . 

A different type of two-component system, standing in between the neutral and non-neutral Coulomb systems, 
are bilayer containing spatially separated positive and negative charges which are most easily to realize in semi- 
conductors by means of doping (electron-hole bilayers). These systems are of high interest because the strength 
of the correlations can be tuned by varying the layer separation d. The interplay of intra-layer and inter-layer 
correlations in classical bilayers has been studied in detail for macroscopic, e.g. Iillj and mesoscopic lfT2l[T3]| 
systems. Quantum bilayers have been treated much less, see e.g. lfT4l [T5l [TH [TTl and are much poorer under- 
stood. In particular, most investigations have considered symmetric bilayers, where the hole to electron mass ratio 
M = mh/me equals one. However, the typical mass ratio in semiconductors is on the order of Af = 3 . . . 10, 
and even exotic materials exist where M reaches 40 ifTSl or even higher values. 

For this reason, in this paper we concentrate on the effect of the mass ratio on crystal formation in quantum 
electron-hole bilayers. Varying M from 1 to 100 at low temperature and high density, we can tune the hole 
behavior from delocalized (quantum) to localized (quasi-classical) while the electrons remain delocalized all the 
time. As was recently observed for bulk semiconductors Q, holes undergo a phase transition to a crystalUne 
state if the mass ratio exceeds a critical value of Alcr ~ 80. Here, we extend this analysis to bilayers where 
Mcr depends on d and the in-layer particle density. To reduce the complexity of the problem, here we will keep 
d fixed. The complicated overlap of correlation and quantum effects of both, electrons and holes, is fully taken 
care of by performing first-principle path integral simulations. We present results for two types of e-h bilayers: a 
mesocsopic system of = 36 particles in a parabolic trap and for a macroscopic system of the same density. 
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2 Model and Parameters 

The physical realization of the mass asymmetric bilayers considered here can be a system of two coupled quan- 
tum wells filled with electrons and holes, respectively. An additional in-plane potential can produce the lateral 
confinement of the carriers leading to a system of two coupled quantum dots. Recently, we have analyzed in de- 
tail a possible realization of a parabolic in-plane potential using the idea of the quantum Stark confinement |fT9l . 
An inhomogenous electric field applied perpendicular to the QW plane changes the energy of a particle in the 
quantum well because the penetration of a particle inside the barrier material depends on the strength of the elec- 
tric field. For example, in GaAs and ZnSe based QW one can achieve harmonic trap frequencies from 1 GHz to 
1 THz for typical electric field strengths of 10 — 20 kV/cm. 

In this paper, we approximate two coupled QWs by a model of two vertically separated 2D layers popu- 
lated with Ne electrons and Nh holes (we consider the case ~ Nh = N/2). The charges interact via the 
Coulomb potential. The underlying Hamiltonian is well defined and is of practical importance for semiconductor 
heterostructures 

where the electrons (e) are confined to the plane 2 = and the holes (h) to the plane z — d; also and rj are 
the in-plane 2D radius vectors describing the particle coordinates in each layer. In the following all lenghts will 
be given in units of the effective Bohr radius gb = ?i^e/m*e^. For example, for GaAs and ZnSe quantum wells 
this results in the length units aB{GaAs) = 9.98 nm and aB{ZnSe) = 3.07 nm. Energies and temperatures are 
measured in Hartree units: lHa[GaAs) = 11.477716^(133.1 K) and lHa{ZnSe) = 53.937r7ey (625.8 K). 

For the mesoscopic trapped system the density is controlled by the harmonic trap frequency (we use rn*^Lo1 = 
m'^ojf/) and is characterized by the coupling parameter A — (e^/e^o) / (huje) ~ Iq/o-b with /q = h/mloje- In this 
case, the coupling parameter for the holes is related to the electron coupling as A/j = A(777^/77i*)^/''. Also for 
Coulomb systems in a parabolic trap one can find the following usefull relations. For two classical particles in 
a parabolic trap their separation distance tq in the ground state is given by: e^/ero = meijjlr^/^.. Now we can 
define the density parameter (in analogy to the Brueckner parameter ~ {r)/aB for macroscopic systems) 
as follows: 7=^ = ra/aB = {26^ / emeUilY^^ / aB = 2^/'^A'*/'^. We will use this formula to obtain approximate 
relations between the densities in the mesoscopic and macroscopic system by relating A <-> ^ . 

2.1 Numerical details 

To solve the problem of N interacting particles described by the Hamiltonian ([TJ we use the path integral Monte 
Carlo (PIMC) method. The applied PIMC simulation technique was described in detail in Ref. |20|. The effective 
interaction potentials used in the expressions for the high-temperature pair density matrices were obtained by 
using the matrix squaring technique 12111221 . 

One of the main obstacles that limit applicability of the PIMC method for systems of particles obeing Fermi 
statistics is the so called Fermion sign problem. Without additional approximations the direct fermionic PIMC 
simulations are only limited to problems where the degeneracy is not very high. This, certainly depends on the 
physical situation and is related to the particle density, interaction strength and temperature. Full inclusion of 
the quantum exchange effects for the number of particles considered here, i.e N^i^i^^ « 36 — 64, will not be 
possible without neglecting the spin statistics and permutations in the electron and hole subsystems. However, 
direct comparison of the PIMC simulations without spin |23| with the results of Ref. |25| which include spin 
effects show, that the errors introduced by neglecting the spin statistics are of the order of few percents and are 
completely negligible for A > 10. The considered here electron densities, i.e « 18, are sufficiently low (for 
the holes the corresponding parameter 7-i''^ is even larger due to their larger mass) and the dominant effect for 
the interparticle correlations (and in particular for the holes) are driven mainly by the strength of the Coulomb 
interaction and not by quantum statistics effects. Also, there is no doubt that the spin will have a negligible 
effect on the localized states of the holes when they form a 2D lattice. Hence, we expect, that the solid-liquid 
transition investigated in this paper will not be sensitive to the particle spin. Nevertheless, the question about the 
true ground state of the electron liquid (i.e. spin polarized or unpolarized), just after the solid-liquid transition is 
currently under active discussion 11261 and requires further investigation. 
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In the simulations presented below we assume that the electrons and holes can reach thermal equilibrium and 
are cooled down by using, e.g. ^He/'^He dilution refrigeration to a temperature of fc^T = 1/3000 Ha. For 
ZnSe (GaAs) this corresponds to an absolute value of T = 208.6(44.37) mK. At these low temperatures thermal 
fluctuations are negligible and the system is practically in the ground state. 

In the PIMC representation of the density matrix applied in our simulations ||20l . we have used 256 (in some 
cases 128) beads (high temperature factors). This was sufficient to reache convergence for the full energy better 
than 1% and an even better accuracy for the pair distribution functions. The use of such a moderate number 
of beads was only possible by using pre-computed tables of the pair density matrices for all types of Coulomb 
interactions, i.e. for the intra-layer and inter-layer interaction terms in the hamiltonian ([T]i, and for the external 
parabolic confinement for both electrons and holes. To reduce the enormous computational effort for a simulation 
of fermions, here we used Boltzmann statistics for both electrons and holes, and the spin effects are omitted. For 
reasons discussed above we expect that this will not influence the results of this paper significantly. 

Both layers are treated as pure 2D layers of zero thickness. Considering that the thickness of real physical 
QWs is of the order of few Bohr radii, this approximation seems to be reasonable for the range of densities 
considered here, i.e. = {r) /as > 10, and an inter-layer distance of d = 20aB- [The case when d and 
(r) become comparable to the well width would require essentially more computationally costly 3D simulations 
and inclusion additional terms related to the QW potential in the hamiltonian ([T])]. For quite narrow QWs with a 
thickness of about lag and less, the addiabatic approximation can be succesfully used with the 2D hamiltonian ([T]| 
with slightly changed interaction terms (see Ref. Il27l ). 

For the chosen inter-layer distance, d — 2QaB, our system represents essentiall a 3D structure, as the intra- 
layer and inter-layer correlations are on the same length and energy scales. For small ratios d/r^ <C 1 the 
system approaches the single layer limit, whereas for d/rg ^ 1, it behaves like two uncoupled layers. Also, for 
d = 20aB we can completely neglect the inter-layer tunneling, and the life time of electrons and holes can reach 
a few microseconds which is much larger than their equilibration time. 

The mass ratio of the electrons and holes is varied in the range of 1 < m^/ rn* < 100 which covers practically 
all semiconductor materials. In our simulations we have found that the initial equiUbration time needed to bring 
the system from an initial randomly chosen configuration to the thermodynamic one, depends on the electron- 
hole mass ratio and the strength of the external confinement. Usually we skip the first 10 000 — 100 000 MC-steps 
and only then start to accumulate thermodynamic averages. 



2.2 Calculated quantities 

2.2. 1 Pair and radial distribution functions 

The physically relevant quantities to investigate a phase transition are the radial, n{r), and pair distribution 
function, g{r). Both functions are a good probe of the short and long-range order in the system and yield 
information on the importance of correlation effects. In Statistical Mechanics these quantities are given by the 
expressions 

^''^''W^^Tl^EEWI'-'^- 1 n,(,)(r) = -^ 5](<5(|r,|-ro)), (2) 
" i=i j=i i=i 

where a and b are two particle species, Tq is the reference point for the radial density (e.g the center of the 
parabolic potential as used here), and (. . .) denotes the thermodynamic average. In the PIMC approach the 
averaging is performed with the TV— particle density matrix, i.e 



dridr2 ■ . . dr at (. . .) p(ri,r2, . . . ,rAr;/3). (3) 



After the high-temperature decomposition this integral includes also additional integrations over the particle 
coordinates on the intermediate "time-slices" and, as a result, the particle images on each time slice also contribute 
to the distribution function which significantly improves the convergence of the simulations. 
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2.2.2 Lindemann parameter 

One of the criteria to investigate structural phase transitions (e.g., solid-liquid phase transition) was proposed 
by Lindemann |28|, who used vibration of atoms in the crystal to explain the melting transition. The average 
amplitude of thermal vibrations increases with temperature of the solid. At some point the amplitude of the 
vibrations becomes so large that the atoms start to occupy the space of their nearest neighbors and disturb them, 
and the melting process is initiated. According to Lindemann, the melting might be expected when the root mean 
vibration amplitude \J {6v?) /a^ exceeds a certain threshold value {{Su'^) is the particle fluctuation from a lattice 
site, a — l/y/Tm, n is the density). Namely, when the amplitude reaches at least 10% of the nearest neighbor 
distance, this quantity exhibits a rapid growth when the temperature becomes close to the melting temperature 
of the solid phase. While for 3D systems this criterion can be successfully used, in 2D this quantity shows a 
logarithmic divergence, \n{L/a), with the increase of the system size L. Instead, to indicate the phase transiton 
from a liquid to a crystal, in 2D, one should apply the modified Lindemann criterion and use the relative distance 
fluctuations L29J 

where ry is the distance between the particles i and j. To reduce the effect of particle diffusion through the 
cluster (in a finite system) or through the simulation cell (for a macroscopic system), which leads to very slow 
convergence with the increase of the system size, in the calculation of ^ we have performed partial averaging 
over 1 000 MC-steps (one block). After the current block has been completed we proceed to a new one and the 
MC averaging was repeated for the next 1 000 MC-steps. The difference in the fluctuations measured from block 
to block can characterize the ordering in the system and is more effective for large systems. 



2.2.3 Nature of the phase transition in 2D systems 

Strictly speaking, in classical macroscopic 2D systems at T 7^ a true crystal state does not exist. The absence of 
off-diagonal long range order in the system manifests itself in the existance of two disordered phases characterized 
by different asymptotic behavours of the pair correlation function g{r, r'). The system undergoes a transition at 
a finite temperature Tkt (Kosterlitz-Thouless transition) when the asymptotes g{r, ?'')||r_r'|^oo changes from 

( _ exp(-|r- r'|/^(T)) X ^ ( 1 (rr^rr ^ ... 

ai^^"-)- |r„rr(^) iT>TKT) to g(r,r)«^— -^p^ (T < T^t). (5) 

The important question of the relevence of the standard Kosterlitz-Thouless theory also for 2D quantum systems 
has been disscussed in Ref. II 3 21 for the two-dimensional XY model. A generalization for Coulomb systems is 
subject of ongoing work ll33ll . Concerning the interpretation of the results of the present publication we indeed 
find a abrupt transition in the decay of the maxima and mimima of g(r, r') (see the discussion below) which can 
be approximated by the asymptotes in Eq. (|5]l. 



3 Numerical results: Mesoscopic system 

In the following we consider a bilayer system populated with a mesoscopic number of = Nh — 36 electrons 
and holes. The results of our simulations are presented in Figs.|2ll3] In our simulations two different densities are 
analyzed, given by A = 5 and A = 10.5 which corresponds to the first maxima of the pair distribution function 
= 8.7 and gjj'jf'^ = 19, respectively. These densities are chosen such that, at the given temperature, hole 
crystallization is expected to occure, at least for large mass ratios M = mh/rrie- If the density is chosen too low, 
the Coulomb coupling would to weak for crystallization. On the other hand, if the density is too high, the crystal 
vanishes due to quantum melting. At the chosen densities the electrons are always in the quantum liquid-like 
state, while the state of the holes can be changed by varying M. 

At the chosen densities the total cluster radius is i?A=5 = TOas (i?A=io.5 = ISOas). That means that the 
average densities (in a single layer) are for GaAs 77a=5 = 9.4 • lO^/cm^ and 71a=io.5 = 2.0 • 10^/cm^, and for 
ZnSe nA=5 = 9.9 • 10^"/ cm^ and ^^=10. 5 = 2.2 • lO^'^/cm^. These values are for the electrons, for the holes the 
radius slightly decreases when M is increased. 
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Fig. 1 (Color online) Path integral Monte Carlo configuration of holes (blue points) and electrons (red dots) in a bilayer 
system with distance d = 20as, temperature T = //a/3000 and different mass ratios: M = 5 (left column), M = 20 
(center) and M = 100 (right). Each particle is represented by 256 dots (path integral) which, for the electrons, are mutually 
penetrating. First two rows: 36 electrons and holes in a harmonic trap with coupling strengths A = 5 (upper panel) and 
A = 10.5 (second panel). Shown is a typical snapshot (without statistical averaging). Note the different axis scales in the 
two panels. Row 3: Simulation snapshots of a macroscopic bilayer (Ne — Nh = 64 electrons and holes in the simulation 
cell with periodic boundary conditions, the borders mark the simulation cell. Each particle is shown only once). The density 
matches the one in the confined system of the second row. There are structural defects as the triangular lettice is not uniform. 



Consider first Fig. [T] which gives an overview on the observed behavior for the two densities (first two rows) 
when the mass ratio is varied in the range from 1 to 100. The first observation is that, in all cases, the electrons 
are distributed almost continuously, whereas the holes become localized when M exceeds 20 (5) at A 5 (10.5). 
Due to the rotational symmetry of the trap, the holes are arranged in concentric shells. 

The main difference between the mesoscopic system with a parabolic inplane confinement and an infinite 
system are well-known finite size effects, see e.g. [34| which are related to the rotational symmetry instead of 
translational symmetry. Further, even when averaged over the modulation caused by the shells the density is not 
constant over the entire system, cf. left part of Fig. |2l The average density is highest in the center and decreases 
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radial in plane distance R [a„] hole-hole distance r [a ] 



Fig. 2 (Color online) Hole radial distribution (left fig.) and hole-hole pair distribution (right) for A = 5 and five mass ratios 
(see inset) for a mesoscopic confined bilayer with A^'e = Nh = 36. 

towards the cluster surface. Fig. |2] also clearly shows the effect of the mass ratio. With increasing M the hole- 
hole correlations increase leading to increased hole localiczation |5 6|. This is accompanied by a pronounced 
modulation of the radial density n{R) and the pair distribution (PDF) 17/1^, see Fig.|2] The reduction of the zero 
point fluctuation with increase of the particle mass M leads to a hole localization and crystal formation. It is 
found that the shell radii in the radial density profile n{R) in Fig. |2] as well as the peak positons in the hole- 
hole pair correlation function ghh of the mesoscopic cluster (A — 5) are independent from mass ratio M. For 
M = 100 we find that the holes are arranged in 3 shells populated with 16, 12, 7 and a single particle in the 
center, see Fig.[T] 

Fig. |2] shows, that by changing the mass ratio from 1 to 100 the holes exhibit a transiton form a delocalized 
quantum state with wave function overlap to a highly ordered quasi classical state, while the electrons stay in 
a quantum fluid state and their correlations change only little with M for the present parameters. We note that 
the classical Coulomb coupling parameter for M ~ 100 is Fa^s = {Ucorr) / {Ukin) = 345 and rA=io.5 — 158, 
which is beyond the critical value for the macroscopic (OCP) crystallization Tcrit = 137. 

Let us now consider the response of the electrons to the formation of the hole crystal. While the electron 
density is almost structurelss, some details can be seen in the electron-hole PDF, Fig. |3] This function has a 
distinct peak at zero (in-plane) distance showing the electrons and hole are pairwise vertically aligned for afl 
values of M. Also, the next peaks of the e-e PDF are aligned with those of ghh- The small shift in the peaks of 
the two functions is due to the normalization. In order to compare the details of the cluster arrangements with 
the macroscopic system below, in Fig. [3] we have divided g^h and g^h by the corresponding functions sor M = 1 
where they are structureless. This allows to largely eliminate the effect of the trap (but slightly shifts the extrema). 

4 Numerical results: Macroscopic system 

To understand the relevance of our above mesoscopic results for larger systems contaning hundreds or thousands 
of particles we performed additional simulations for a macroscopic e-h bilayer without confinement potential. We 
have considered — Nh ~ 64 electrons and holes in a simulation cell of the size {L^ x Ly} = {76.185aB x 
65. 97803} with periodic boundary conditions (PBC). This corresponds to a density parameter « 10 (average 
particle distance in units of the electron Bohr radius). This density was chosen to be comparable to the average 
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Fig. 3 Hole-hole (upper fig.) and electron-hole (lower fig.) pair distribution functions for the mesoscopic bilayer with 
Ne = Nh = 36 and A = 10.5 for four values of M (see inset). The curves are normalized to the corresponding PDF for 
the case M = 1 to eliminate the influence of the decay of the average density in the trapped system, cf. Fig. |2] Note that 
electrons and holes are always pairwise aligned vertically. The plot includes distances up to two times the radius which causes 
the increase of the PDF for large distances. 



density in the finite system (see Sec.O for the case of couphng parameter A = 10.5. The mass ratio M was 
varied between 1 and 100, the temperature was fixed to kT = l/3000i?a. The number of particles and the 
dimensions of the cell, i.e Ly = ^/iLx/l, were choosen to best fit the symmetry of a triangular lattice, which is 
expected to be formed by the holes. We note that finite size effects are of the order of few percents, a systematic 
analysis with larger particle numbers is beyond the scope of this paper. 

Let us now consider the results for the macroscopic bilayer. Three typical shapshots for M — 5, 20, 100 
are shown in the lower row of Fig. [T] As in the mesoscopic system, for all cases the electrons are completely 
delocalized. In contrast, the hole localization increases from i\/ = 5 to M = 100. Also, we confirm that the 
density of the mesoscopic system (second row) is well matched: the average distance between two holes as well 
as their extension (given by the size of the blue dots) is very close to the trapped case. 

Consider now the pair distributions. In Fig.|4](upper fig.) we show the hole-hole PDF for different mass ratios 
1 < Af < 20. Since the particle number and box size is fixed, the average particle density stays constant and 
the position of the first peak of the PDF are practically independent of M. However, the general behavior of the 
PDF changes drastically. For A/ > 4 we observe clear oscillations typical for the solid phase. Even the third and 
fourth peaks are well resolved (the scale exceeds half of our simulation box). These oscillations become rapidly 
damped by changing Af to 3 and below, here the PDF show liquid-like features. The third peak is now strongly 
suppressed. This transition can be quantified by computing the ratio of the (magnitude of the) first minimum to 
the first maximum which is 71 = 0.48, for M = 4, and 71 = 0.65, for M — 3. Similarly, for the third peak this 
ratio becomes 73 = 0.76 and 73 = 0.96, respectively. The ratio 71 is frequently used as an empirical criterion for 
the solid-liquid transition in classical systems; in a one-component 3D system the critical value is known to be 
7]^ « 1/3. If a universal values exists also in the present two-component 2D quantum system where the transition 
is expected to be of the Kosterlitz-Thouless type is an interesting question which deserves further analysis. 
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inplane distance [a^] 

Fig. 4 Hole-hole (upper fig.) and electron-hole (lower fig.) pair distribution functions for the macroscopic bilayer with 
Ne = A'^h = 64 (with periodic boundary conditions) for the mass ratios M = 1, 3, 4, 5, 10, 20 (the maxima increase with 
increasing M). Note the alternating location of maxima and minima of g^h and geh. 



Let US now compare the pair distributions with those in the mesoscpic system at the same density (A = 10.5), 
Figs. |4]and|3] Interestingly, we find that the first peaks of g^h have approximately the same height, and also 
the peak positons are very close, see upper parts of the two figures. Further we observe that the minima of guh 
are significantly deeper in the macroscopic case. This is explained by intershell rotations which occur in the 
mesoscopic system |5| and wash out the correlations. The present results are at temperatures above the freezout 
of these rotations. 

Consider now the relative importance of the inter-layer correlations for the stability of the hole crystal. To this 
end, we have plotted the e-h PDF in Fig.|4](lower fig.). For the symmetric case, M = 1, and also for Af = 3 we 
observe similar behavior: the highest probability has the configuration where the electrons reside (in their own 
layer) just below the holes, as was observed in the mesoscopic system, lower part of Fig. |3] Obviously, the height 
of this peak is small, the modulation depth is around 2% because of the high electron degeneracy (delocalization). 
This means that these peaks cannot be associated with bound states (indirect excitons) since the electron density is 
well above the Mott density riMott for this system where excitons break up because the repulsion of two excitons 
exceeds the electron-hole binding. Note that riMott depends on the layer separation d which governs the binding 
energy and the typical size a^^ of an indirect exciton which is of the order of d. Hence, for the present parameters, 
d/aB — 20 and w 10, the in-plane exciton size exceeds the separation of two neighboring electrons which 
causes exciton ionization. On the other hand, reducing d below 10, excitons become stable (for temperatures 
below the exciton binding energy) which is confirmed by our PIMC simulations. 

For larger mass ratios, M > 4, a completely different behavior of g^h emerges. From our analysis of the 
hole-hole correlations we know that the holes are now in an "ordered state" (or, in the terms of the Kosterlitz- 
Thouless theory, in a "less disordered state" with a power-low decay of off-diagonal long-range order). Now, 
there is no maximum of g^h at zero distance, and the function exhibits oscillations. The explanation is that 
the electron density is modulated due to the presence of the hole crystal with maxima located in between the 
holes. While the amplitude of the oscillations is small, (about 1% modulation depth) they are clearly visible 
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Fig. 5 (Color online) Left Fig.: Relative distance fluctuations of the holes, u^^ , (Eq. as a function of the mass ratio M 
for a macroscopic (black solid line) and mesoscopic bilayer system for two densities (see inset). Right Fig.: Decay of the 
amplitude of the maxima and minima of the hole correlations {\ghh — 1|) in the macroscopic system, cf. Fig.|4l for the seven 
mass ratios - from bottom to top: Ad = 1,2,3, 4, 5, 10, 20. Note the change from an exponential (for AI < 3) to a power law 
decay (for M > 4) which signals the Kosterlitz-Thouless transition. 



and become systematically more pronounced when M increases, see Fig. |4] (lower part). We, therefore, expect 
that appearance (disappearance) of these oscillations of geh is an additional indicator of a phase transition in the 
present asymmetric bilayer system. 

Finally, as another quantity sensitive to phase transitions, we consider the relative distance fluctuations u!^^ 
of the holes, Eq.lH as a function of M, Fig. |5] (right part). This quantity exhibits a rapid drop between M = 3 
and M — A which is related to a localization transition. We can translate from the critical mass ratio (which is 
expected to be between three and four) to the hole density parameter ri'^\ using ri''' = r'f^ mh / and the 
position of the first peak of the hole-hole PDF at Ts ~ 10. As a result, we obtain that the phase transition in 
the hole layer occurs at a critical density in the range 30 < Ts < 40. This result is close to the value r* « 37 
known as the critical density of solid-liquid transition in the one-component quantum 2D system at T = ll24l . 
Compared to this value, in our bilayer system, we observe indications of stabilization of the "ordered state" of 
the holes due to presence of the electron layer. 

We note that, at smaller values of d (e.g. d = bas and d — lOas) no hole crystal is found. Instead we 
observe formation of indirect excitons which form a solid phase of composite particles. At the same time, the 
interparticle interaction changes from Coulomb to dipole-like which reduces the value of the classical coupling 

2 .2 

parameter to F = |^//cbT. Similar tendencies have also been also in simulations of symmetric classical and 
quantum e-h-bilayers lfm[T4ll . 



5 Discussion 

Analyzing the peak height (amplitude) of the g^h in the macroscopic system (Fig. |4|l in dependence on the peak 
positon rfih we can deduce to the correlation decay law and compare to the asymptotics (|5]l. In the disordered 
phase of small mass ratios M = 1...3 we find an exponential correlation decay of ghh, see right part of Fig.|5] 
From mass ratio M = 1 to Af = 3 disordering is lowered and the correlation length increases from ^ = 6.5 to 
^ = 13. Increasing the mass ratio above the critical mass ratio, i.e. M > 4, we find a topological transition to the 
Kosterlitz-Thouless phase with power law correlation fall-off. 

We may now obtain a critical mass ratio at which quantum melting of the hole crystal takes place. Using as a 
criterion a critical value of = 0.15 of the relative hole-hole distance fluctuations we obtain M^riti^ = 5) ~ 5 
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and Mcrit(A — 10.5) « 2.8, in the mesoscopic system, and Mcrit{rs = 10) ~ 3.1, in the macroscopic system 
(recall that it corresponds to A w 10.5). Obviously, the absolute numbers are somewhat arbitrary, but the allow 
for an analysis of the dominant trends, i) Merit depends on density. It decreases when the coupling strength A 
increases in agreement with earlier observations for small e-h clusters ||231 . ii), there is good agreement between 
the critical mass ratios of the mesoscopic and the macroscopic system (within 10%). iii) the critical values are 
much smaller than the value of Merit ~ 80 in a 3D bulk system [7J which underlines the remarkable additional 
control of physical behaviors existing in a bilayer system by a variation of the layer separation d. It is expected 
that further reduction of d will allow to further reduce Merit and to increase the maximum density of the hole 
crystal to values below r^s^^ = 20 1131 . 
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